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ABSTRACT 

We have measured the weak gravitational lensing of faint, distant background galaxies by MS 1054-03, a rich 
and X-ray luminous cluster of galaxies at a redshift of z = 0.83, using a two-colour mosaic of deep WFPC2 images. 
The small corrections for the size of the PSF and the high number density of background galaxies obtained in these 
observations result in an accurate and well calibrated measurement of the lensing induced distortion. 

The strength of the lensing signal depends on the redshift distribution of the background galaxies. We used 
photometric redshift distributions from the Northern and Southern Hubble Deep Fields to relate the lensing signal 
to the mass. The predicted variations of the signal as a function of apparent source magnitude and colour agrees 
well with the observed lensing signal. The uncertainty in the redshift distribution results in a 10% systematic 
uncertainty in the mass measurement. 

We determine a mass of (1.2 ±0.2) x IQ'^^ hl^ M© within an aperture of radius 1 /igy Mpc. Under the assumption 
of an isothermal mass distribution, the corresponding velocity dispersion is 131l;!;g|. For the mass-to-light ratio 
we find 269 ± 37/!5oM(^Lb0 after correcting for pass-band and luminosity evolution. The errors in the mass and 
mass-to-light ratio include the contribution from the random intrinsic ellipticities of the source galaxies, but not 
the (systematic) error due to the uncertainty in the redshift distribution. However, the estimates for the mass and 
mass-to-light ratio of MS 1054-03 agree well with other estimators, suggesting that the mass calibration works 
well. 

The reconstruction of the projected mass surface density shows a complex mass distribution, consistent with 
the light distribution. The results indicate that MS 1054-03 is a young system. The timescale for relaxation is 
estimated to be at least 1 Gyr 

We have also studied the masses of the cluster galaxies, by averaging the tangential shear around the cluster 
galaxies. Using the Faber-Jackson scaling relation, we find the velocity dispersion of an L* galaxy (L^ = 8 x 
1O1°/iJ^Lb0 for MS 1054-03) is 203 ± 33 km/s. 

Subject headings: cosmology: observations - gravitational lensing - dark matter - galaxies: clusters: individual 
(MS 1054-03) 



1. INTRODUCTION 

Studying the mass distribution of massive clusters at high 
redshifts is important as their existence puts strong constraints 
on possible cosmological models (e.g. Eke et al. 1996; Bah- 
call & Fan 1998; Donahue et al. 1998). Furthermore, these 
systems may be young and thus much can be learned about the 
formation of massive clusters of galaxies. 

Massive structures in the universe induce a systematic dis- 
tortion of the images of faint background sources, providing a 
powerful diagnostic to study their mass distributions (e.g. Mel- 
lier 1999). The weak lensing signal allows us to map the pro- 
jected mass surface density of the lens (e.g. Kaiser & Squires 
1993), and determine its mass. Although the weak lensing 
mass estimate does not require assumptions about the dynam- 
ical state or geometry, it does require a good estimate for the 
redshift distribution of the faint sources. 

Since Tyson, Wenk, & Valdes (1990) succesfully measured 
the systematic distortion of faint galaxies by a foreground clus- 
ter, many massive clusters of galaxies have been studied (e.g. 



Bonnet, Mellier, & Fort 1994; Fahlman et al. 1994; Squires et 
al. 1996). These studies concentrated on massive clusters with 
redshifts z- 0.2-0.5. 

Small et al. (1994) tried to measure the weak lensing sig- 
nal for the optically selected, high redshift cluster CI 1603-1-43 
(z = 0.89). Its observed velocity dispersion is 1226^[5| km/s 
(Postman, Lubin, & Oke 1998). Small et al. (1994) used source 
galaxies down to a limiting magnitude of / = 25, but failed to 
detect a signal. Their results are consistent with a lower mass 
for the cluster. This is supported by the cluster's X-ray lumi- 
nosity of L,(0. 1-2.4 keV) = 1 . 1 x Hy^^-l ergs/s (Castander et 
al. 1994), which is an order of magnitude lower than the X-ray 
luminosity of the brightest, high redshift. X-ray selected clus- 
ters 

The X-ray emission of clusters of galaxies provides an ef- 
ficient means to select massive clusters of galaxies. Since the 
work of Small et al. (1994), several high redshift. X-ray lu- 
minous clusters have been discovered from Einstein (Gioia & 
Luppino 1994) and ROSAT observations (Henry et al. 1997; 
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Ebeling et al. 1999; Rosati et al. 1998). 

Luppino & Kaiser (1997; LK97 hereafter) successfully de- 
tected the weak lensing signal from the X-ray selected, high 
redshift cluster MS 1054-03 using deep ground based imag- 
ing. Clowe et al. (1998) reported on the detection of weak 
lensing signals for the X-ray selected, high redshift clusters 
MS 1137H-66 and RXJ1717H-67, from deep Keck imaging. 
However, these studies suffered from the lack of knowledge of 
the redshift distribution of the faint sources. 

Lensing studies of high redshift clusters are difficult because 
the lensing signal is low and most of the signal comes from 
small, faint galaxies. This is the regime where HST observa- 
tions have great advantage over ground based imaging, because 
the background galaxies are much better resolved. As a result, 
higher number densities of sources can be reached, the images 
are not crowded, and more importantly the correction for the 
circularization by the PSF is much smaller Especially for high 
redshift clusters, where the redshift of the lens approaches that 
of the sources, the lensing signal depends strongly on the red- 
shift distribution of the sources. 

In this paper we present our weak lensing analysis of the 
cluster MS 1054-03 using a mosaic of deep WFPC2 images. 
It is the second cluster we studied using a mosaic of deep HST 
observations. In Hoekstra et al. (1998, HFKS98 hereafter) we 
presented the analysis of CI 1358H-62, a rich cluster of galaxies 
at a redshift of z = 0.33. 

MS 1054-03 (z = 0.83) is the most distant cluster in the Ein- 
stein Extended Medium Sensitivity Survey (Gioia et al. 1990; 
Gioia & Luppino 1994). It is extremely rich and has a restframe 
X-ray luminosity' Lv(2- 10 keV) = 2.2 x lO'^^/jjo ergs/s (Don- 
ahue et al. 1998), which makes it one of the brightest X-ray 
clusters known. From their ASCA observations Donahue et al 
(1998) find a high X-ray temperature of 12.3!^;^ keV. 

Our space based observations allow us to study the cluster 
mass distribution in much more detail than LK97. Photomet- 
ric redshift distributions inferred from the Hubble Deep Field 
North (Fernandez-Soto, Lanzetta, & Yahil 1999) and South 
(Chen et al. 1998) give a good approximation of the redshift 
distribution of the sources used in this weak lensing analysis. 
As a result it is possible to obtain well calibrated mass esti- 
mates for high redshift clusters of galaxies from multi-colour 
imaging, even if the availability of spectroscopic redshift data 
for the sources is limited. 

The observations and data reduction are outlined in section 2. 
In section 3 we briefly discuss the method we use to analyse the 
galaxies. Section 4 deals with the light distribution of the clus- 
ter. In section 5 we present the measurements of the lensing 
induced distortions of the faint galaxies and the reconstruction 
of the projected mass surface density. Section 6 deals with the 
redshift distribution of the faint galaxies. Our mass determina- 
tion is presented in section 7. In this section we also discuss 
several effects that complicate accurate mass determinations. 
The results for the mass-to-light ratio are given in section 8. In 
section 9 we examine in more detail the substructure of the clus- 
ter. The lensing signal due to the cluster galaxies themselves is 
studied in section 10. 

2. OBSERVATIONS 

In this analysis we use WFPC2 images taken with the Hub- 
ble Space Telescope (GO proposal 7372, PI Franx). Figure [l| 
shows the layout of the mosaic constructed from the 6 pointings 

'Throughout this paper we will use ^50 = Hf/{50 knv'^Ts/lpc), 0^=0.3 and Ha = 



of the telescope. The integration time per pointing was 6400s 
in both the F606W and F814W filter. For the weak lensing 
analysis we omitted the data of the Planetary Camera because 
of the brighter isophotal limit. The total area covered by the 
observations is approximately 26 arcmin^. 
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Fig. 1 . — Layout of the observed field. The six pointings are 
indicated and numbered. The positions of the 20 brightest clus- 
ter galaxies are indicated by the black dots. The dashed circle 
indicates an aperture with a radius of 1 /i^q Mpc. The image is 
oriented such that North is up. The area covered by the obser- 
vations is approximately 26 arcmin"^. 



For each pointing the exposures of a given passband were 
split into two sets. The images in each set were offset from 
each other by integer pixels (±3 pixels). Both sets were re- 
duced and combined separately, resulting in two images (for 
details see van Dokkum et al. 1999; van Dokkum 1999). These 
two images were designed to be offset by 5.5 pixels, to allow 
the construction of an interlaced image (e.g. Fruchter & Hook 
1998; van Dokkum 1999), which has a sampling that is a factor 
V2 better than the original WFPC2 images. 

Except for the images in the first pointing in F814W, all off- 
sets were within 0.1 pixel of the requested value. We found that 
we had to omit the shape measurements in the F814W band 
for pointing 1 , because of the deviating offset for this exposure. 
However, we could use the F6Q6W image, and as the noise in 
the shear estimate is dominated by the intrinsic ellipticities of 
the galaxies, the loss in signal-to-noise is minimal. 



3. OBJECT ANALYSIS 

Our analysis technique is based on that developed by 
Kaiser, Squires, & Broadhurst (1995; KSB95 hereafter) and 
LK97, with a number of modifications which are described in 
HFKS98. 

The first step in the analysis is the detection of the galaxy im- 
ages, for which we used the hierarchical peak finding algorithm 
from KSB95. We selected objects with a significance i/ larger 
than 5(7 over the local sky. The detected objects were analysed, 

0. This gives a scale of 1" = 9.5 Kl kpc at the distance of MS 1054. 
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yielding estimates for the sizes, magnitudes and shapes of the 
objects. The detection and analysis were done per WFPC2 chip 
and filter As a new feature, our software computes an estimate 
for the error on the polarization from the data. This error esti- 
mate is discussed in Appendix A. 

The resulting catalogs were inspected visually in order to re- 
move spurious detections, like diffraction spikes, HII regions 
in resolved galaxies, etc. 4402 objects were detected in the 
F6Q6W images, whereas the F814W images yielded 4290 stars 
and galaxies. The combined catalog consisted of 5848 stars and 
galaxies, of which 2844 were detected in both filters. 
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For all detected objects we estimated the apparent magnitude 
using an aperture with a diameter of 6v2'"g pixels (where rg is 
the Gaussian scalelength of the object as given by the peak find- 
ing algorithm). We converted the measured counts to F6Q6W 
and F814W magnitudes, zero-pointed to Vega, using the zero- 
points given in the HST Data Handbook (STScI, Baltimore). 
We added 0.05 magnitude to the zero-points to account for the 
Charge Transfer Efficiency (CTE) effect. 

Figure g shows the apparent magnitude versus the half light 
radius of the object. Moderately bright stars are easily identi- 
fied as they are located in the vertical locus. Bright stars satu- 
rate and hence their measured sizes increase. Using figure g we 
selected a sample of moderately bright stars, which were used 
to examine the PSF. 

To obtain a catalog of galaxies we removed the stars from 
the catalog. At faint levels the discrimination between stars and 
galaxies is not possible, but by extrapolating the number counts 
of moderately bright stars to faint magnitudes, we find that the 
stars contribute less than 2% to the total counts. 

We attempted to determine colours for all detected objects 
(including galaxies detected in only one filter), using the same 
aperture in both filters. The software failed to determine colours 
for some of the galaxies. Most of these were extremely faint, 
and detected only in the F814W band. For the weak lensing 
analysis we only used galaxies for which colours were deter- 
mined. Due to this selection we lost 30% of the galaxies fainter 
thanF814W = 26.5 (533 galaxies). 

The shape parameters of the galaxies were corrected for the 
PSF (cf. section 3.1) and the camera distortion^ We used the 
estimated errors on the shape measurements (cf. Appendix A) 
to combine the results from the F6Q6W and FE14-W images in 
an optimal way, resulting in a final catalog of 4971 galaxies. 

The better sampling of the interlaced images improved the 
object detection and analysis significantly. HFKS98 showed 
that due to the poor sampling of WFPC2 images only the shapes 
of galaxies with sizes rg > 0"12 could be measured reliably. 
The images used here allow us to use all detected galaxies. 
HFKS98 found a number density of 79 galaxies arcmin"^ from 
3600s exposures in both F606W and F8 14W. The observations 
presented here yield a number density of 191 galaxies arcmin"^. 
This higher number density is due to both the longer integration 
time, and the improved sampling of the interlaced images. An 
analysis of the effect of small pointing errors in the interlaced 
images shows that these are not important for our analysis (cf. 
Appendix B). 

Figure ^ shows the observed number counts of galaxies in 
our observations versus apparent magnitude. For comparison 
we also show the counts from the HDF-North (Williams et al. 
1996). This indicates that our catalog of galaxies is complete 
down to F6Q6W = 26.5 and F814W = 26. 



3.1. Corrections 

Several observational effects systematically distort the im- 
ages of the galaxies. To obtain an accurate and unbiased es- 
timate of the weak lensing signal these effects have to be cor- 
rected for To do so, we follow the scheme outlined in KSB95, 
LK97, and HFKS98. 

^In HFKS98 we discussed the correction for tlie sliear introduced by the camera. It was subsequently pointed out by Rhodes, Refregier, & Groth (1999) that the 
camera shear in fact is radial, and not tangential as presented in HFKS98. Because the effect is small, and we used a mosaic of images, it does not change the results 
of HFKS98. 



Fig. 2. — Plot of the apparent magnitude in the F606W filter 
(upper panel) and the F814W filter (lower panel) versus half 
light radius. In both figures the vertical locus where stars are 
found is easily identified. 
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Fig. 3. — Plot of the number counts of detected galaxies in 
the F6Q6W band (solid circles) and F814W band (open trian- 
gles). The dashed line shows the number counts from the HDF- 
N (Williams et al. 1996) in the F814W filter, and the soHd fine 
shows the counts in the F606W filter. This indicates that our 
catalog is complete down to F6Q6W = 26.5 and FE14W = 26. 
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Fig. 4. — The polarization field of stars taken from the in- 
terlaced observations in the F814W filter. The orientation of 
the sticks indicates the direction of the major axis of the PSF, 
whereas the length is proportional to the size of the anisotropy. 
The polarization values are measured using a Gaussian weight 
function with a dispersion of 0''07. The lower left panel corre- 
sponds to chip 2, lower right to chip 3, and the upper right one 
denotes chip 4. We have omitted chip 1, which is the planetary 
camera. 



HFKS98 showed that the PSF of WFPC2 is highly 
anisotropic at the edges of the chips. We selected a sample 
of moderately bright stars from the observations of MS 1054. 
Figure ^ shows the observed polarization of these stars in the 
F814W band. The anisotropy pattern is similar to the pattern 
presented in HFKS98. 

Because the PSF changes slightly with time, we fitted the 
shape parameters of our new observations to a scaled model 
of the PSF polarization measured from the globular cluster M4 
and added a first order polynomial: 

pr = a-pf + co + cix + C2y. 

This procedure gave excellent fits to the observed PSF 
anisotropy in our observations of MS 1054. 

The PSF also circularizes the images of the small, faint 
galaxies, thus lowering the amplitude of the lensing signal. 
Compared to ground based observations, the small faint galax- 
ies are much better resolved in HST observations, resulting in 
much smaller corrections for the circularization by the PSF. 
This allows us to measure a well calibrated lensing signal. 

To obtain estimates for the real shapes of the faint galaxies, 
we determined the 'pre-seeing' shear polarizability P'' (LK97; 
HFKS98). The measurements of P''' for individual galaxies are 
rather noisy, and therefore we binned the data as a function of 
size of the object. 
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Fig. 5. — Required correction factor VP'^ to estimate the 'pre- 
seeing' ellipticity of the object. The solid line shows the re- 
sults for the F606W filter, and the dashed line shows the cor- 
rection factor for F814W. The pre-seeing shear polarizability 
P^ corrects for the circularization by the weight function and 
the PSF. It depends strongly on the size of the object, and in- 
creases rapidly as the object's size approaches the size of the 
PSF. As the PSF in ground based data is much larger, the cor- 
rections for galaxies of similar sizes would be much larger. The 
upper panel shows the cumulative number of objects as a func- 
tion of size. The objects for which the peak finder assigned a 
size smaller than the PSF are analysed with a fixed weight func- 
tion matched to the PSF. These objects end up in the bin of the 
smallest objects. 
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The results are presented in figure ^ which shows the re- 
quired correction factor to obtain the true elHpticity of a galaxy 
from the observed polarization. The upper panel shows the cu- 
mulative number of objects as a function of object size. The ob- 
jects for which the peak finder assigned a size smaller than the 
PSF were analysed with a weight function matched to the PSF. 
These objects end up in the bin corresponding to the smallest 
objects. 

The correction for objects with sizes comparable to the size 
of the PSF is large. The PSF in ground based images is much 
larger, and therefore large correction factors are required for ob- 
jects that are still well resolved in HST images. From a 0."65 
deep ground based image we found that P'^ w 0.2 for galaxies 
of R ^ 25, whereas galaxies of similar brightness in our HST 
images have P'' w 0.6. 

As a test of our procedure we compared the distortion mea- 
sured from the F6Q6W images and the F814W images (which 
included the correction for the circularization) using only ob- 
jects for which shape parameters could be measured in both 
images. We compared the average tangential distortions and 
found that the results for the two filters agree well to one an- 
other (we found {gj) = 0.0366 ± 0.0066 from the F6Q6W data 
and (gr) = 0.0375 ± 0.0068 for the F814W images). 

Recently it has been found that the Charge Transfer Effi- 
ciency problem for WFPC2 has worsened. This can system- 
atically change the shapes of the faint galaxies, as charge is 
'smeared' along the CCD columns. This could result in a de- 
tectable change in the shapes of the galaxy images. 

We compared the shapes of the objects in the regions where 



the pointings overlap (cf. figure [|). Unfortunately the num- 
ber of objects was low, and the contribution of shot noise to 
the shape estimates substantial. We did not detect a systematic 
change in the shapes, but we cannot put strong constraints on 
the size of the effect. We also simulated the effect, assuming 
that the change in shape increased linearly with row number 
for each chip, with no change at low row numbers and a shear 
of 10% at high row numbers. We found that this resulted in a 
neglegible change in our mass estimate. 

4. LIGHT DISTRIBUTION 

Figure ^ shows a plot of the colour of the galaxies versus 
their F814W magnitude. The cluster is visible as the concen- 
tration of galaxies with F6Q6W - F^lAW ~ 2. The galaxies on 
this colour-magnitude relation are generally early type galaxies, 
which are among the reddest objects in the cluster. 

We selected 324 galaxies brighter than F814W = 24.5, lying 
close to the colour-magnitude relation, to study the cluster light 
distribution. This catalog was compared to the spectroscopic 
catalog (van Dokkum 1999). The contamination by field galax- 
ies is small, and we removed galaxies with redshifts incompati- 
ble to that of the cluster. Many of these belong to a small group 
at z = 0.547 (van Dokkum 1999). 

We have use SExtractor (Bertin & Arnouts 1996) to deter- 
mine total magnitudes for these galaxies. Figure^ shows grey 
scale plots of the luminosity and number density weighted light 
distributions. The images have been smoothed with a Gaussian 
with a FWHM of 20 arcseconds. The galaxy number density 
and the luminosity weighted distributions look quite similar. 
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Fig. 6. — Colour-magnitude diagram of the objects for which colours could be determined. The cluster colour-magnitude relation is 
clearly visible. The errorbars at the bottom of the plot correspond to the errorbars on the F814W magnitude. 
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Fig. 7. — (a) Luminosity distribution of cluster galaxies brighter than F814W 24.5. The contours are 10%, 20%, etc. of the peak 
value; (b) Number density distribution of this sample of cluster galaxies; (c) Number density distribution of all detected galaxies 
brighter than F814W 25; (d) Number density distribution of all detected galaxies fainter than F814W 25, where we corrected the 
counts for the area covered by the (large) bright galaxies. In panels b-d the interval between subsequent contours are 20 galaxies 
arcmin"^, and the thick lines indicate levels of 100 and 200 galaxies arcmin"^. All distributions have been smoothed with a Gaussian 
filter with a FWHM of 20 arcseconds (indicated by the shaded circle). The mosaic covers approximately 26 arcmin"^, and the total 
field shown is 7 by 6 arcmin. 
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The cluster light distribution is elongated east-west, and three 
clumps can be identified. 

Figure ^ also shows gray scale images of the smoothed num- 
ber density when we split the sample of all detected galaxies 
into a bright (F814W < 25) and a faint sample (F814W > 25). 
The cluster is clearly visible in the bright sample, but does not 
show up in the number counts of the faint galaxies (which have 
been corrected for the area covered by large bright galaxies). 
Figure ^ shows several significant peaks, indicative of some 
clustering at high redshift. 

We estimate the luminosity in the redshifted B band follow- 
ing van Dokkum (1999). The direct transformation from the 
HST filters to the passband corrected B band is given by 

B, = FSUW-0.03iF606W-FSUW)+1.23, 

where B^ denotes the passband corrected B band magnitude. 
The luminosity is given by 

Lb = IqO'^('^bq-B:+DM+Afshw)t ^ 

where Mbq = 5.48 is the solar absolute B magnitude, DM is 
the distance modulus, and Afg,i4w is the extinction correction 
in the FEUW filter towards MS 1054-03. The redshift of 0.83 
for MS 1054 gives adistance modulus of 44. 09-51og/i5o. Tak- 
ing galactic extinctions from Burstein & Heiles (1982) we use 
a value of 0.03 forAf^i4w- 

At faint magnitudes, the colour-magnitude relation is less 
clearly defined. To estimate the total light of the cluster, an 
estimate of the amount of light contributed by galaxies fainter 
than F814W = 24.5 is required. To do so, we fitted a luminosity 
function to the sample of bright, colour selected cluster galax- 
ies. We found that a Schechter function with a = -1 .0 and M^ = 
-21.8 + 51og/!5o (L* = 8 X 10"'/iJ^Lbq) could fit the observed 
counts of the sample of bright, colour selected cluster galaxies. 
For nearby clusters one typically finds M^ = -21.0 + 51og/!5o 
(eg. Binggeli, Sandage, & Tammann 1988), indicating that the 
galaxies in MS 1054-03 are approximately twice as bright as 
nearby cluster galaxies. This is in good agreement with the re- 
sults of van Dokkum et al. (1998) who studied the fundamental 
plane in MS 1054-03. The luminosity function indicates that 
we miss approximately 5% of the light when not including the 
faint galaxies. 

MS 1054-03 is a high redshift cluster and one might expect 
a relatively high fraction of blue cluster galaxies, which are 
too blue to be included in our colour selected sample of clus- 
ter galaxies. Comparison of our sample to a sample of bright, 
spectroscopically confirmed members (van Dokkum 1999) in- 
dicates that we miss 16% of the total light in the B band. 

The cumulative light profile as a function of radius from the 
centre for the sample of bright galaxies is presented in figure H. 
The profile is multiplied by a factor 1.05 x 1.16 to account for 
the light from the faint cluster galaxies and the bluest galaxies. 
The solid line corresponds to the observed profile, whereas the 
dashed line shows a isothermal profile for comparison. We es- 
timate the total luminosity within an aperture of 1 /i^q Mpc to 
be 1.0 X lO^^h-^lLBQ. 



5. WEAK LENSING SIGNAL 

Due to the intrinsic ellipticities of the galaxies, their individ- 
ual shape estimates are very noisy measurements of the lensing 
induced distortion. We therefore need to average the measure- 



ments of a large number of galaxies to obtain a useful estimate 
for the distortion g. 

To estimate the average distortion, we weight each object 
with the inverse square of the uncertainty in the distortion, 
which includes the contribution of the intrinsic ellipticies of 
the galaxies and the shot noise in the shape measurement (cf. 
Appendix A): 
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Fig. 8. — The cumulative, passband corrected, B band luminos- 
ity as a function of radius from the assumed cluster centre. The 
solid line corresponds to the observed profile. The profile was 
calculated using the colour selected sample of cluster galaxies, 
and was corrected for the contribution of very blue galaxies and 
faint galaxies. The dashed line corresponds to an isothermal 
profile, which is shown for comparison. The arrow indicates a 



radius of 1 /150 Mpc. 
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Selecting galaxies with 22 < F814W < 26.5 and ^606^^- 
F814W < 1.6 (the 'source' sample from table |l]), and smooth- 
ing using a Gaussian with a FWHM of 20 arcsec, we obtain 
the distortion field presented in figure ^. The sticks indicate 
the direction of the distortion, and the length is proportional to 
the amplitude of the distortion. It shows that on average the 
faint galaxies tend to align tangentially to the (elongated) clus- 
ter mass distribution, as expected from gravitational lensing. 

The average tangential distortion gj in radial bins is a useful 
measure of the lensing signal. It is defined as gj = -{g\ cos2(/)+ 
g2sin2(/)), where (j) is the azimuthal angle with respect to the 
assumed centre of the mass distribution. We assume that the 
position of the Brightest Cluster Galaxy defines the centre of 
the cluster At large radii from the cluster centre, where we 
do not have data on complete circles, we account for possbile 
azimuthal variation of gj by fitting 

gr(r»= (§r)W(l+§2,cCos(20)+g2,.,sin(2(/.)) (2) 

to the observed distortions. This procedure gives the correct 
value for an ellipsoidal mass distribution. Although the cluster 
appears elongated, no large azimuthal variations in the tangen- 
tial distortion were found. 
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Fig. 9. — Distortion field g obtained from the sample of source 
galaxies, smoothing the data using a Gaussian with a FWHM 
of 20 arcsec (indicated by the shaded circle). The orientation 
of the sticks indicates the direction of the distortion, whereas 
the length is proportional to the amplitude of the signal. Due 
to gravitational lensing the images of the faint sources tend to 
align tangentially to the cluster mass distribution. 



The average tan gen tial distortion as a function of radius is 
presented in figure p^ . In the absence of lensing, the points 
should scatter around zero. The figure shows a systematic tan- 
gential alignment of the sources. The dependence of the signal 
on radius is quite complicated, indicative of a complex mass 
distribution. The results of the two-dimensional mass recon- 
struction (cf. section 5.1) show that the sharp decrease in the 
tangential distortion around 40" is due to substructure in the 
cluster mass distribution. 

As a test, figure p^ shows the results when the phase of the 
distortion is increased by 7i/2 (i.e. rotating the sources by 7i/4). 
In this case the signal should vanish if it is due gravitational 
lensing, as is observed. 

We divide our sample of galaxies in several subsamples, 
based on their apparent magnitude and colours. Table |l] lists 
some details of the subsamples. It shows that even our 'bright' 
sample in fact consists of rather faint galaxies. 

The tangential distortion profiles of the various subsamples 
are presented in figure O. A clear lensing signal is detected 
for all subsamples. The signals of the bright and the faint sam- 
ple are independent measurements, and agree well with one- 
another Although the profiles of the blue and the red sample 
agree well, the difference in the amplitude of the lensing sig- 
nal is clear This is due to the fact that, on average, the redder 
galaxies are at lower redshifts than the blue galaxies. In sec- 
tion 6 the issue of the redshift distribution of the background 
galaxies is discussed in detail. 
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Fig. 10. — (a) Tangential distortion as a function of radius from the cluster centre using the sample of source galaxies as defined in 
table [l|. A clear lensing signal is detected, although the dependence of the signal on radius indicates a complex mass distribution, (b) 
Measured signal when the phase of the distortion is increased by 7i/2. If the signal observed in (a) is due to gravitational lensing, gx 
should vanish as is observed. In both figures the arrows indicate a radius of 1 /jjq Mpc. 
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Fig. 1 1 . — Tangential distortion as a function of radius from the 
cluster centre for the various subsamples: (a) bright galaxies; 
(b) faint galaxies; (c) blue galaxies; (d) red galaxies. The sig- 
nals for the bright and the faint galaxies are quite similar in both 
strength and shape. Although the shapes of the signals for the 
blue and the red sample are similar, the difference in strength of 
the signal is clear. In these figures the arrows indicates a radius 
of 1 /jg,) Mpc. 



5.1. Projected mass distribution 

Figure |l2| shows the reconstructed mass surface density 
based on the distortion field presented in figure ^. We used 
the maximum probability extension of the original KS algo- 
rithm (Kaiser & Squires 1993; Squires & Kaiser 1996), which 
does not suffer from some of the biases of the original KS93 in- 
version, and can be applied directly to fields with complicated 
boundaries. Several effects complicate a direct interpretation of 
the mass map (cf. section 7.1). However, if the distortions are 
relatively small, it gives a fair description of the cluster mass 
distribution. 

The mass reconstruction shows that the cluster mass distri- 
bution is quite complex: three distinct peaks are visible. The 
positions of the peaks coincide fairly well with the peaks in the 
light distribution (cf. Fig. ^. The results presented in LK97 
lacked the resolution to detect the substructure. 



Donahue et al. (1998) studied the X-ray properties of 
MS 1054 and found that the distribution of the gas is elongated 
similar to the weak lensing mass distribution. They also find 
evidence for substructure. However, the peak in their map of 
the X-ray brightness coincides with the western clump. Deeper 
X-ray observations by Neumann et al. (1999) show that the 
peak in the X-ray is close to the central clump. In section 9 we 
investigate the substructure in more detail. 

Towards the edges of the observed field the noise in the re- 
construction increases rapidly, and therefore the significance of 
the structures at the edge of the field is low. 

The reconstruction after increasing the phase of the distor- 



tion by 71/2 is shown in figure 123. As this corresponds to the 



imaginary part of the mass surface density, no signal should be 
present if the observed distortion is due to gravitational lensing. 
We find no significant large scale structures in this map. 

6. REDSfflFT DISTRIBUTION 

To relate the observed lensing signal to a measurement of the 
mass one has to estimate the mean critical surface density 

Scrit = c'(47rGD,/3)-', 
where D/ is the angular diameter distance to the lens and 

/3=(max(0,D,/D,)), 

where Dj and D;, are the angular diameter distances from ob- 
server to the source and lens to the source respectively. The an- 
gular diameter distance to MS 1054-03 is 1.96 /iJq Gpc, which 
yields Scit = 845/?"'/i5oM0pc"- (using 0„, = 0.3, and ^Ia = 0). 

The question that remains is the value of (/?). The uncer- 
tainty in the redshift distribution of the sources has hampered 
the mass estimates from previous lensing studies of high red- 
shift clusters of galaxies (LK97; Clowe et al. 1998). 

Down to / '^ 23 the redshift distribution of galaxies is fairly 
well known from redshift surveys, but the galaxies used in 
weak lensing analyses are generally too faint to be included in 
spectroscopic surveys. However, photometric redshift distribu- 
tions derived from broad band photometry of the Northern HDF 
(Fernandez-Soto, Lanzetta, & Yahil 1999) and Southern HDF 
(Chen et al. 1998) may be used to obtain an estimate of {f3) as 
a function of apparent magnitude and colour. The photometric 
redshift distributions can be compared to our data, because we 
use the same filters, and the catalogs are complete down to the 
magnitude limit of our weak lensing analysis (F814W < 26.5). 
Figure 13 shows the comparison of the weak lensing signal with 
the predicted trends based on the HDFs in independent magni- 
tude and colour bins. 
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(3) 


(4) 


(5) 


(6) 


name 


F814W 


F606W-F814W 


# galaxies 


n 


m 




[mag] 


[mag] 




arcmin"^ 


[mag] 


source 


22-26.5 


<1.6 


2643 


102 


25.2 


bright 


22-25 


<1.6 


1137 


44 


24.2 


faint 


25-26.5 


— 


1906 


73 


25.8 


blue 


22-26.5 


<0.8 


1198 


46 


25.4 


red 


22-26.5 


0.8-1.6 


1445 


55 


25.1 



Table 1: Properties of the various subsamples taken from the final catalog of galaxies. 
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The value of (/3) does not only depend on the redshifts of the 
sources, but it also depends on the cosmological parameters that 
define the angular diameter distances. Figure p^shows how the 
average value of f3 depends on il,„ and Oa for MS 1054. For 
a given value of the cosmological constant, the dependence on 
ftm is weak, but given i7„,, the average value for f3 is more sen- 
sitive to the value of the cosmological constant. In an J7a domi- 
nated universe the value of (3 is approximately 20% higher than 
in a low f^A universe, and as a result the estimated mass for 
MS 1054-03 would be 20% lower 

In our estimates of the lensing signal we weight the contribu- 
tion of each source according to its error on the distortion (cf. 
appendix A). Thus faint galaxies, which tend to be at higher 
redshifts are given less weight. Applying the same weight func- 
tion (as a function of apparent magnitude) to the HDF galaxies 
shows that the effective average (3 is slightly (6%) lower than 
the value indicated in figure |l^. We use (/3) = 0.23 (flA = 0) for 
our catalog of sources. Using 17,,, = 0.3 and ft^ = 0.7 yields a 
higher value of 0.27 for (/3) . 

7. MASS ESTIMATES 

One of the advantages of gravitational lensing over other 
methods to estimate masses of clusters is that no assumptions 
have to made about the dynamical state or geometry. From the 
weak lensing signal we can directly estimate the projected mass 
within some aperture, although this requires rather good esti- 
mates for the redshift distribution of the background sources 
(cf. section 6). 

However, several effects complicate a precise measurement 
of the total mass of a cluster. These are discussed in section 7.1. 
In section 7.2 we present our mass estimate for MS 1054-03. 
The result is compared to other estimators in section 7.3. 



Fig. 12. — (a) Reconstruction of the projected mass surface 
density. Towards the edges of the observations the signal-to- 
noise in the reconstruction decreases rapidly. The origin co- 
incides with the assumed cluster centre. The reconstruction 
shows that the mass distribution consists of three distinct mass 
concentrations, similar to the light distribution, (b) The recon- 
structed mass surface density when the phase of the distortion 
is increased by 7i/2. As this projects out the imaginary part of 
the surface density, the signal should vanish if it is due to grav- 
itational lensing. The mass distribution was smoothed using a 
Gaussian with a FWHM of 20" (indicated by the shaded circle). 
The interval between adjacent contours is 0.05 in k. 

Though there is some variation between the redshift distri- 
butions in the Northern and Southern fields, the results from 
the weak lensing analysis agree well with the predictions, giv- 
ing us confidence that the calibration of the lensing signal using 
photometric redshifts works. 

The average photometric redshift ofF^lAW < 26.5 galaxies 
in the HDF South is higher than it is in the Northern field. If 
this difference is representative of field to field variations the 
redshift distribution the uncertainty in weak lensing mass esti- 
mates of high redshift clusters (z ^ 0.8) is significant: on the 
order of 10% (a bootstrapping analysis of the photometric red- 
shift catalog indicates that the relative uncertainty in (/3) for 
each individual deep field is only 2%). For clusters at lower 
redshifts the introduced uncertainty is smaller. 



7.1. Complications in accurate mass determinations 

7.1.1. Mass sheet degeneracy 

The mass sheet degeneracy (Gorenstein, Shapiro, & Falco 
1988) reflects the fact that the observed distortion is invariant 
under the transformation: 



k' = (\-\)k+\. 



(3) 



Thus scaling the surface density, while adding a sheet of 
constant surface density, leaves the (observed) distortion un- 
changed. Measuring the distortion out to large radii, and argu- 
ing that the cluster surface density is neglegible at large radii, 
solves the problem if no intervening sheet of matter is present. 
In the analysis presented here we assume that this is the case. 

The degeneracy can be lifted by measuring the magnification 
of the sources relative to an unlensed population of sources (e.g. 
Broadhurst, Taylor, & Peacock 1995). However, it is difficult 
to measure the magnification effect with sufficient accuracy to 
detect a mass sheet A of a few percent. 

7.1.2. Measuring large distortions 

Simulations of sheared images show that the correction 
scheme we use to measure the shapes of the galaxies works 
well for relatively small distortions (< 0.3). The KSB95 ap- 
proach assumes that the distortion is small everywhere. The 
pre-seeing shear polarizability P'^ is derived from the average 
over all objects, and hence it is not accurate in the regime where 
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Fig. 13. — (a) Comparison of the photometric redshift distribu- 
tion with the observed weak lensing signal. The solid line is the 
predicted value of (/3) as a function of apparent magnitude in 
the F814W band, using the photometric redshift distribution of 
galaxies with F606W-F814W< 1.6 and F814W < 26.5. The 
dashed line shows the results using only the photometric red- 
shifts from the HDF-North, and the dotted line gives the results 
for the Southern field. The variation in the value of (/3) obtained 
from the two deep fields introduces an additional uncertainty in 
the weak lensing mass estimate. The points correspond to the 
weak lensing signal, which is obtained by a SIS model fit to the 
data (K(r) = r£/2r, where r^ corresponds to the Einstein radius.) 
to the data at radii larger than 75 arcseconds. The right y-axis 
indicates the vertical scale in rg- (b) The comparison as a func- 
tion of colour of the background galaxies. We used galaxies 
with 22 < F814W < 26.5. For the reddest bins (open circles), 
the contamination by cluster members is significant (30-50%), 
resulting in an underestimate of the lensing signal. The points 
in both figures are independent, and clearly show that the red 
galaxies carry a lower lensing signal and thus their average red- 
shift is lower than that of the blue galaxies. The dependence 
of the lensing signal on the apparent magnitude is quite weak 
for the faint galaxies. In both figures the lower panels show a 
histogram of the number of objects in each bin. The widths of 
the bins are indicated by the horizontal errorbars. 




1.0 



Fig. 14. — Contour plot of (/3) as a function of 17,,, and J^a = 1, 
based on the photometric redshift distributions inferred from 
the HDF-N (Fernandez-Soto et al. 1999) and HDF-S (Chen et 
al. 1998). For a given value of 51a, the value of (/3) depends 
only slightly on ft,,,. 

the distortions are large. When the average distortion is 0.4, the 
method already underestimates the distortion by 10%, and the 
discrepancy increases with increasing distortion. Other meth- 
ods, like that proposed by Kuijken (1999), are more accurate 
when the distortion is large. 

Due to the high redshift of MS 1054-03 the distortions are 
relatively small (cf. figure nOb, and the KSB95 still works. 
However, placing the cluster at z = 0.2-0.4 would produce 
much larger distortions, and the effect mentioned here has to 
be corrected for. 

7.1.3. Converting g to j 

The data provide a measure of the distortion g, which is re- 
lated to the shear by 7 = (1 - K)g. In the inner regions of the 
cluster, where k is significant, the factor (1 - k) is important. 

7.1.4. Broad redshift distribution 

Another effect is related to the redshift distribution of the 
sources. When computing the average distortion, we use an 
average value for /? for the sources, i.e. we assume that the red- 
shift distribution can be approximated by a sheet of sources at 
redshift corresponding to the average (3. As we show now, this 
results in an overestimate of the lensing signal. 

For an individual galaxy the distortion is 



7 



A7o 



1- 



1 - (3sKo 



(4) 



where Koo and 700 correspond to the dimensionless surface den- 
sity and shear for a source at infinite redshift, and Ps = (y(3oo ■ 
If we adopt a common value of f3 for all background galaxies, 
we estimate the distortion as 



1-(/3.)kc 



(5) 



12 



Weak lensing analysis of MS 1054-03 



whereas the true average distortion is in fact given by the aver- 
age of equation Q: 



is)'' 



A7o 



I-/3.5K0 



(6) 



The ratio g^g indicates by what factor the shear 7 is overesti- 
mated. To first order it can be approximated by: 



(A 
(1) 



1+ 



f (/?.?) 



V(A> 



1 K=l + 



(^ 

(/3)^ 



1 ]k. 



(7) 



This result was also derived by Seitz & Schneider (1997) in 
their analysis of non-linear mass reconstructions. Using the 
photometric redshift distributions for the HDFs, we find that 
for MS 1054-03, one would overestimate the shear by a factor 
(1 + 0.62k): this represents a 7% effect at a radius of 500 /15Q 
kpc. The size of the effect increases with lens redshift as is 
demonstrated in figure Ol 




Fig. 15. — (a) Plot of {{(5^)/{PY- 1) as a function of lens red- 
shift, based on the distribution of photometric redshifts from 
the HDF-N and HDF-S. To produce this figure we used galax- 
ies with F6Q6W-FUm/ < 1.6 andF814W < 26.5. The effect 
increases with lens redshift. (b) Plot of {{j3'^)/{j3Y- 1)k as a 
function of lens redshift. We used the value of k at a distance 
of 0.5/150 Mpc for a cluster with an isothermal mass distribu- 
tion and a velocity dispersion of 1000 km/s. For MS 1054-03 
the shear at this distance from the cluster centre would be over- 
estimated by 7%. 



7.1.5. Magnification 

A further source of error arises from the fact that lensed im- 
ages of background sources are not only distorted, but also 
magnified. Their total flux is increased by a factor /i = ((1 - 
k)^ - 7^)"' . As a result the average redshift of the sources in- 
creases. The change in j3 is approximately l.lKdfVdm, which 
is small. 

All mass estimates presented below have been corrected for 
the effects listed above. Since they are proportional to k, they 
mainly affect the inner regions of the cluster. 



7.2. Mass of MS 1054-03 

To obtain an estimate of the mass of MS 1054, we fit a sim- 
ple mass model to the data. It is straightforward to include 
corrections for the effects mentioned in section 7.1.1 into the 
model. However, the result depends on the assumed radial sur- 
face density profile, and given the complex mass distribution of 
MS 1054, fitting a simple model might not yield the best mass 
estimate. Nevertheless, it should give a fair indication of the 
mass (within a factor of 1.5, depending on the radial profile). 
We fitted a singular isothermal sphere model (k = r£/2r) to the 
tangential distortion at radii larger than 75". We found a value 
of r£ = 1 1 ."5 ± 1 "4, which corresponds to a velocity dispersion 
cr= 1365 ±84 km/s. 

Another approach is often referred to as aperture mass densit- 
ometry. The shear can be related directly to a density contrast. 
We will use the statistic of Clowe et al. (1998), which is related 
to the C statistic of Fahlman et al. (1994): 



Ccin) = 2 f din r{jT)+ 2(1-^)-' f d\nr{jT)- (8) 

<J max •/ 



Cc yields the mean dimensionless surface density interior to 
ri relative to the mean surface density in the annulus from r2 to 

^max- 

Cc(ri) = R(r' < ri)-R(r2 < r' < r^^x)- (9) 

It shows that one can measure the average surface density 
within ri up to a constant if the shear were observable. We used 
the two-dimensional mass reconstruction, to obtain an estimate 
of the shear from the observed distortion. 
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Fig. 16. — Plot of the Cr statistic as a function of radius from 
the assumed cluster centre. The values for (c were calculated 
using the sample of source galaxies (cf. table H), and using 
r2 = 105.2" (1 h'l Mpc) and r^ax = 200". At smafl radii it is 
difficult to correct for the effects mentioned in 7.1.1, and as a 
result these points are overestimates of the true values. The 
arrow indicates a radius of 1 /igy Mpc. 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


n,„ 


^A 


1 Aji Mpc 


(/3> 


Sen, 


M(< 1/i^iMpc) 


O" 






["] 




[r'feoMepc-^] 


[IO'^Aj^Mq] 


[km/s] 


0.3 


0.0 


105.2 


0.23 


845 


1.25±0.17 


mi-li 


0.3 


0.7 


93.8 


0.27 


754 


1.07±0.12 


12151? 


1.0 


0.0 


120.6 


0.23 


969 


1.38±0.19 


1379!^? 



Table 2: Mass estimates for various choices of fl,,, (1) and fl\ (2); (3) angular scale of 1 /i^g Mpc in arcseconds; (4) average /3 for 
our catalog; (5) critical surface density; (6) mass within an aperture of 1 /ijg Mpc radius; (7) corresponding velocity dispersion if the 
mass distribution were isothermal. 



Figure [16| presents the measurement of Cir) using the sam- 
ple of source galaxies as defined in table |I] (cf. Fig. |lO|), and 
taking r^ax = 200". At large radii from the cluster centre we 
cannot average the data on complete circles, and therefore we 
used equation || to estimate the average tangential shear The 
results show that the surface density in the centre is quite high, 
making it difficult to accurately correct the measurements at 
small radii for the effects listed in section 7.1.1. 

We find C- = 0.070± 0.014 at a radius of 1 /ij^ Mpc (indicated 
by the arrow). This yields a lower limit on the mass within this 
aperture of (8. 0± 1.7) x 10'"^ /z^g M©. To obtain the total mass, 
we need to estimate the average surface density in the outer 
annulus. To this end we use the results of the SIS model fit. 
This gives 0.038 ±0.005 for the average dimensionless surface 
density in the outer annulus. Thus we find that the total mass 
within an aperture of 1 /zjg Mpc is (1.2±0.2) x 10^^ /jjg M©. If 
the mass distribution were isothermal, this mass estimate would 
correspond to a velocity dispersion of 1311^^^ km/s. 

Table g lists the mass estimates for different cosmologies. It 
shows that the mass estimate for MS 1054-03 decreases for an 
increasing value of the cosmological constant. 

The uncertainties in the mass estimates presented in this sec- 
tion only reflect the contribution from the intrinsic ellipticities 
of the source galaxies. As shown in section 6, the uncertainty 
in the redshift distribution of the sources may contribute an ad- 
ditional 10% systematic uncertainty. 



7.3. Comparison with other estimates 

For some time MS 1054-03 was the only very massive. X-ray 
selected, high redshift cluster known, and therefore it has been 
studied extensively. 

LK97 were the first to detect a weak lensing signal for this 
cluster. The tangential distortion presented in LK97 is about 
a factor 1 .5 higher than the signal from our HST observations 
(using galaxies with 21.5 < F814W < 25.5). The difference in 
the signals might be due to the large corrections for the circular- 
ization by the PSF in the ground based observations of LK97. 
Furthermore, their mass estimate was hampered by the lack of 
knowledge of the redshift distribution of the faint galaxies. 

Donahue et al. (1998) measured a high X-ray temperature of 
12.3;!;2'2 ^^^ from ASCA observations, which corresponds to a 
velocity dispersion of 1400 ± 170 km/s (90% confidence limits) 
under the assumption of an isothermal mass distribution. This 
estimate is in excellent agreement with our weak lensing mass 
estimate. However, the agreement is somewhat surprising: the 
cluster is clearly not relaxed, and one might expect that the gas 
is shock heated as the three massive clumps are merging. 

The dynamical velocity dispersion was measured by van 



Dokkum ( 1 999), who finds a value of 1 1 50 ± 90 km/s based on 
81 galaxies. Earlier reported measurements based on smaller 
samples (Donahue et al. 1998; Tran et al. 1999) are consistent 
with this result. 

All the evidence gathered so far indicates that MS 1054 is 
a very massive cluster, and using a Press-Schechter analysis 
one can estimate the likelihood of finding such a cluster in the 
EMSS catalog. Based on such analyses, LK97, and Donahue 
et al. (1998) have argued that the presence of MS 1054 in the 
EMSS catalog is hard to reconcile with il„, = 1. However, in an 
open or f^A dominated universe, the likelihood of finding mas- 
sive high redshift clusters is much greater (e.g. Bahcall & Fan 
1998). If the primordial density fluctuations are non-gaussian, 
the presence of MS 1054-03 in the EMSS catalog does not nec- 
essarily exclude a high density universe (Willick 1999). 




40 60 80 

radius ["] 

Fig. 17. — Average mass-to-light ratio within an aperture for 
radii larger than 50 arcseconds from the cluster centre (filled 
circles). The mass in the aperture was estimated using equa- 
tion H and assuming that the average dimensionless surface den- 
sity in the outer annulus k(105'.'2 <r< 200") = 0.038 ±0.005. 
The points are correlated and the errorbars only reflect the un- 
certainty in the measurement of the distortion. The open circles 
denote the observed mass-to-light ratios of the three clumps in 
the mass distribution, using the results listed in table ^ The 
arrow indicates a radius of 1 /150 Mpc. 
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[1013 hllLBQ] 


[/!5oMg/Lb0] 


[huMq/hsQ] 


0.3 


0.0 


1.0 


124±17 


269 ±37 


0.3 


0.7 


1.2 


90 ±10 


225 ±25 


1.0 


0.0 


0.85 


163 ±22 


299 ±41 



Table 3: Average mass-to-light ratio within an aperture with a 1 /j^q Mpc radius for different choices of the cosmological parameters 
(columns (1) and (2)); (3) Total light within the aperture; (4) observed average mass-to-light ratio in the passband corrected B band; 
(5) mass-to-light ratio in the B band, corrected for luminosity evolution as inferred from the evolution of the fundamental plane for 
the various choices of ilm and JIa (van Dokkum et al. 1998). 



8. MASS-TO-LIGHT RATIO 

Combining measurement of the mass with the estimate of 
the total light in the passband corrected B band within an 
aperture of 1 h^Q Mpc radius (cf. section 4) yields (124 ± 
17)/i5()M0/Lb0. This result has to be corrected for luminos- 
ity evolution in order to compare it to values found for lower 
redshift clusters. Studies of the fundamental plane of clusters 
of galaxies at various redshifts have quantified the luminosity 
evolution (e.g. van Dokkum & Franx 1996; Kelson et al. 1997; 
van Dokkum et al. 1998). The fundamental plane of MS 1054 
was studied by van Dokkum et al. (1998). They found that 
the mass-to-light ratio of early type galaxies in MS 1054-03 is 
~ 50% lower than the present day value. Under the assump- 
tion that the total luminosity of the cluster has changed by the 
same amount, the corrected mass-to-light ratio of MS 1054 is 
269±37/!5oMg/Lbo. 

The resulting mass-to-light ratio depends on the cosmologi- 
cal parameters. In table the observed mass-to-light ratios in 
the passband corrected B band are listed (column (4)). Col- 
umn (5) lists the mass-to-light ratio corrected for luminosity 
evolution as determined from the evolution of the fundamental 
plane for MS 1054-03 for the various cosmologies considered 
here (van Dokkum et al. 1998). 

HFKS98 measured a mass-to-light ratio of (175 ± 
26)/i5()M(7/LB0 for CI 1358H-62, which changes to (231 ± 
34)/25()M(7/Lb0 after correction for luminosity evolution. Carl- 
berg, Yee, & Ellingson (1997) studied a sample of 16 rich clus- 
ters of galaxies, and they find a value of (134±9)/i5oM(r/Lr0 
for the average mass-to-light ratio. Using B-r ^ 1.1 as found 
for nearby early type galaxies (J0rgensen, Franx, & Kjaergaard 
1995) this ti-ansforms to (260 ± 17)/!5oM0^b0 in the B band. 



This comparison indicates that the mass-to-light ratio of 
MS 1054 is similar to what is found from other studies, sug- 
gesting that the range in cluster mass-to-light ratios is rather 
small. LK97 argued that MS 1054-03 might be a relatively dark 
cluster, but our analysis does not confirm their results. 

We also studied the radial dependence of the mass-to-light 
ratio using equation M and an estimate for the average surface 
density in the control annulus. This allows us to estimate the 
average mass-to-light ratio M(< ryLB(< r) as a function of ra- 
dius. To avoid confusion with the clumps in the mass distribu- 
tion, we present the results for radii larger than 50 arcseconds in 
figure nl (solid points). The open circles indicate the mass-to- 
light ratios measured for the three clumps in the mass distribu- 
tion, using the results listed in table H The results indicate that 
the light is more concentrated towards the centre of the cluster 
than the mass. 

9. SUBSTRUCTURE 

The data enable us to study the observed substructure in the 
cluster mass distribution in detail. Figure Il8|a shows the mass 
reconstruction of the central 3.3 by 2.7 arcminute of the HST 
mosaic. The origin coincides with the assumed cluster centre. 
Figure n8p is the overlay of the reconstruction and the F8141V 
image. To allow a direct comparison of the mass reconstruction 
with the light distribution, we also present the overlay of the 
smoothed light distribution in figure Il8|c. This shows that the 
projected mass distribution compares well with the Ught distri- 
bution. 

Next we estimate the masses of the three clumps from the 
observed distortions (see Table 0). We use two measures: 

(/) a simultaneous fit of three isothermal spheres to the full 
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0.184 ±0.040 


0.050 ±0.006 


1.7±0.3 
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100 ±17 
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623 ±76 


0.194 ±0.037 
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2.2 


77 ±13 
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628 ±80 


0.181 ±0.037 


0.050 ±0.006 


1.7±0.3 


1.5 


113±19 



Table 4: Mass estimates for the subclumps. (1) position of the clump; (2) velocity dispersion obtained by simultaneously fitting 3 
singular isothermal spheres to the observed distortion; (3) the observed value of C,c at 250 hl^li kpc; (4) the average dimensionless 
surface density in the control annulus (k(100" <r < 150")) assuming an isothermal mass distribution with a= 131 1 ± 85 km/s ; (5) 
total projected mass within an aperture of 250 HIq kpc radius; (6) total luminosity in an aperture of 250 /iJq kpc radius; (7) average 
mass-to-light ratio in an aperture of 250 h^l kpc. 
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distortion field; 
(//) the Cc-statistic (equation]^). 

A simultaneous fit of three SIS models to the observed distor- 
tion shows that the masses of the three clumps are similar, with 
corresponding velocity dispersions of 600-700 km/s. Compar- 
ison with the estimate of the total cluster mass indicates that 
most of the mass is in the three clumps. For the (c statistic 
we need to estimate the average dimensionless surface density 
in a control annulus around each clump (we use r2 = 100" and 
'"max = 150" as at these radii the effect of the substructure is 
expected to be small). To do so, we assume that the surface 
density profile at large radii is isothermal with a velocity dis- 
persion of (7= 1311±85 km/s as found above. The results are 
listed in Table |. 

The results from the C^-statistic correspond to the total pro- 
jected mass in the aperture, which also includes mass from the 
other clumps. Therefore we consider the results of the fitting 
procedure (column (1)) as the best estimates for the mass. 

However, to estimate the average mass-to-light ratio of the 
clumps within a radius of 250 /ijg kpc we do use the Q statistic, 
as both mass and light are measured in the same aperture. The 
mass-to-light ratios are around 90/!5oM(7/Lbo, somewhat lower 
than the average mass-to-light ratio within an 1 /ijy Mpc radius 
aperture around the cluster centre. 

New X-ray telescopes like Chandra or XMM will have the 
sensitivity to study the X-ray properties of high redshift clus- 
ters in detail. MS 1054-03 will be an interesting target as we 
know how the mass is distributed. Thus high resolution mea- 
surements of the X-ray emission and temperature gradients can 
provide important information about the thermodynamics of 
cluster formation. 



9.1. Relaxation timescale 

The redshift catalog of van Dokkum (1999) enables us to 
study the velocity structure of the cluster. The galaxies more 
that 100" west of the centre show a velocity difference of about 
800 km/s with respect to the average velocity. However, their 
contribution to the total velocity dispersion is small. The three 
clumps are located closer to the centre, and cannot be separated 
in velocity. Therefore the radial motion of the clumps appears 
to be small compared to the velocity dispersion of the cluster 

Both outer clumps in the mass distribution are in projection 
500/150 kpc away from the central clump. If we assume that all 
motion is perpendicular to the line of sight, and using a rela- 
tively slow approach (500 km/s) we estimate the timescale for 
relaxation to be on the order of 1 Gyr Therefore the cluster 
should be relaxed at a redshift of z ^^ 0.65. If the physical sep- 
arations of the clumps are larger, the relaxation time increases 
correspondingly. 






10. GALAXY-GALAXY LENSING 

To Study the mass distribution of MS 1054-03, we concen- 
trated on the lensing signal on relatively large scales. How- 
ever, the cluster galaxies themselves introduce a small signal as 
well. Measuring this galaxy-galaxy lensing signal can be used 
to study the halo properties of field galaxies (e.g. Brainerd, 
Blandford, & Smail 1996; Schneider & Rix 1997, Hudson et 
al. 1998). 



Fig. 18. — (a) Gray scale image of the mass reconstruction of 
the central 3.3 by 2.7 arcminute region of the HST mosaic; (b) 
Overlay of the mass map on the F814W image of MS 1054-03; 
(c) Overlay of the smoothed light map and the F814W image. 
All images have been smoothed using a Gaussian filter with a 
FWHM of 20 arcsec. The interval between adjacent contours 
is 0.025 in n. Comparison of the mass reconstruction with the 
light distribution shows that there is a fair agreement between 
the distributions of mass and light. 
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Weak lensing analysis of MS 1054-03 



A robust method to examine the cluster galaxies is to mea- 
sure the average tangential distortion around an ensemble of 
cluster members. The results are not very sensitive to the pre- 
cise cluster mass distribution, provided one confines the analy- 
sis to small radii. Due to the clustering of the lenses, the signal 
at large radii is dominated by the smooth cluster signal. When 
measuring the signal out to large radii, a careful removal of the 
cluster signal is required. We subtracted the contribution of the 
smooth cluster mass distribution from the observed distortion 
using the two-dimensional mass reconstruction, and computed 
the shear 7 taking the (1 - k) correction into account. 

We used a sample of galaxies brighter than F814 = 24.5 
(~ 0.05L*), with colours close to the cluster colour-magnitude 
relation to study the masses of the cluster galaxies. To mini- 
mize the effect of the large scale cluster mass distribution, we 
excluded galaxies near the centres of the three clumps, result- 
ing in a sample of 241 galaxies. In total, 4300 cluster galaxy- 
background galaxy pairs were used to measure the distortion. 

Massive galaxies produce a larger distortion than less mas- 
sive galaxies and therefore we normalized the amplitudes of 
the signals from the various galaxies correspondingly. To do 
so, we assumed that the shear 7 ex ct^ ex v^Lb, i.e. we assumed 
a Faber-Jackson scaling relation. The results are presented in 
figure |l^. 

A fit of a singular isothermal model to the data, yielded a 
detection of the lensing signal at the 99.8% confidence limit. 
We found for the Einstein radius of an L* galaxy (L^g = 8 x 
10'°/!5gLBQ for MS 1054-03) a value of rg = 0.26 ±0.09 arc- 
seconds, which corresponds to a velocity dispersion of cr = 
203 ± 33 km/s. The mass-to-light ratio within a radius of 20/jJg 
kpc is M/Lb = 7.5 ±2.4/!5o]VKLb0. Correcting this result for 
luminosity evolution yields M/Lb = 17 ± Shs^M/LBo- 

If we neglect the effect of the smooth cluster contribution, 
we find a slightly higher mass for the galaxies: a = 212 ± 32 
km/s. This indicates that the mass estimate determined using 
the direct averaging method is indeed robust. 

Natarayan & Kneib (1998) studied the halo properties of the 
cluster galaxies in AC 1 14. They used a maximum-likelihood 
analysis to constrain the mass and the extent of the galaxy ha- 
los. Constraining the extent of the halos of the cluster galaxies 
in MS 1054 is a delicate task which has not been attempted 
here. It requires a careful analysis of the smooth cluster mass 
distribution in a maximum likelihood analysis. 



an overestimate of the lensing signal. Furthermore, the method 
used (KSB95; LK97) underestimates the distortion when it is 
larger than ~ 0.3. Although we correct for these effects, it is 
best to restrict the analysis to the weak lensing regime, where 
the convergence k is low, and the distortions are small. 
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11. CONCLUSIONS 

We have measured the weak gravitational lensing of faint, 
distant galaxies by MS 1054-03, a rich cluster of galaxies at 
z = 0.83. The data consist of a two-colour mosaic of 6 interlaced 
images taken with the WFPC2 camera on the Hubble Space 
Telescope. 

The expected lensing signal is low for high redshift clusters, 
and most of the signal comes from small, faint galaxies. Com- 
pared to ground based observations, these galaxies are much 
better resolved in HST observations. High number densities 
of sources are reached, and the correction for the circulariza- 
tion by the PSF is much smaller. This enables us to measure 
an accurate and well calibrated estimate of the lensing induced 
distortion of the faint background galaxies. 

Several effects complicate an accurate analysis when the dis- 
tortions or the surface density are non-neglegible (e.g. in the 
cluster cores). We do not measure the shear, but the distortion 
g. Neglecting the effect of a non-zero surface density results in 



Fig. 19. — (a) Average tangential shear around an ensemble 
of cluster galaxies. The signal induced by the smooth cluster 
mass distribution is removed from the observed signal, using 
the mass reconstruction. The signal for each lens galaxy has 
been scaled to match the signal for a L» (L^b = 8 x 1O'°/!5oLb0) 
galaxy, assuming that mass scales with VLb- The solid line is 
a singular isothermal sphere fitted to the measured points. It 
corresponds to a velocity dispersion of cr = 203 ± 33 km/s. (b) 
Average tangential shear around the cluster galaxies when the 
phase of the shear is increased by 7i/2. 

Converting the lensing signal into a mass estimate requires 
knowledge of the redshift distribution of the faint sources. We 
have explored the use of photometric redshift distributions in- 
ferred from the Northern (Fernandez-Soto et al. 1999) and 
Southern (Chen et al. 1998) Hubble Deep Fields. Comparison 
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of the results for the two fields suggest field to field variations 
which lead to a 10% systematic uncertainty in the mass estimate 
for MS 1054-03. Our results for the mass and mass-to-light ra- 
tio of MS 1054-03 agree well with other estimators, suggesting 
that these photometric redshift distributions are a fair approx- 
imation of the true distribution. This gives us confidence that 
we are able to measure a well calibrated weak lensing mass for 
this high redshift cluster 

Using the (^ statistic (Clowe et al. 1998), we find a mass 
(1.2 ±0.2) X lO'^/i5(')M0 within an aperture of 1 h'^^ Mpc ra- 
dius. Under the assumption of an isothermal mass distribution, 
this corresponds to a velocity dispersion of Olliligg km/s. A fit 
of an isothermal sphere model to the observed tangential dis- 
tortion at large radii, yields a velocity dispersion of 1365 ± 84 
km/s. However, this result depends on the assumed radial sur- 
face density profile. The uncertainties in the mass estimates 
only reflect the contribution from the intrinsic ellipticities of 
the sources. The weak lensing mass estimate is in good agree- 
ment with the X-ray estimate from Donahue et al. (1998) and 
the observed velocity dispersion (Tran et al. 1999; van Dokkum 
1999). 

The observed average mass-to-light ratio in the passband 
corrected B band is 124 ± IT/zsoMg/LBo (measured within a 
1 /zjo Mpc radius aperture). We use the results from van 
Dokkum et al. (1998) to correct our estimate for luminos- 
ity evolution. The resulting mass-to-light ratio of 269 ± 
37 /i5oM(7/Lbo can be compared to results obtained for nearby 
clusters, and is in excellent agreement with the findings of Carl- 
berg et al. (1997). 

Our high resolution mass reconstruction shows three dis- 
tinct peaks in the mass distribution, which are in good agree- 



ment with the irregular light distribution. The results of LK97 
lacked the resolution to detect this substructure. In projection 
the clumps are approximately 500 /ig,, kpc apart, and we esti- 
mate the timescale for virialization to be at on the order of ^ 1 
Gyr Thus the cluster should be virialized at a redshift z ^ 0.65. 
We found that the clumps are of roughly equal mass, with cor- 
responding velocity dispersions of ~ 700 km/s. 

To study the masses of the cluster members, we have mea- 
sured the average tangential shear around a sample of 241 
bright cluster galaxies. Assuming the Faber-Jackson scaling 
relation we find a velocity dispersion of 203 ±33 km/s for 
an U galaxy (Ub = 8 x lO'^/iggLB© for MS 1054-03). The 
mass-to-light ratio within a radius of 2Qh'^l^ kpc is M/L^ = 
7.5 ± 2.4/!5oM/Lb0. Correcting this result for luminosity evo- 
lution yields M/Lb = 17 ± 5/i5oM/Lb0. 

The number of known massive high redshift clusters of 
galaxies is increasing rapidly. The results presented here 
demonstrate that these systems can be studied with sufficient 
signal-to-noise to be of interest for studies of cluster forma- 
tion. The improved sampling and throughput of the Advanced 
Camera for Surveys, that will be installed on the HST in the 
near future, will be very useful for weak lensing analyses of 
high redshift clusters of galaxies, especially when the results 
are combined with forthcoming X-ray observations. The com- 
bination of detailed weak lensing analyses and X-ray observa- 
tions by Chandra and XMM will provide important information 
on the thermodynamics of cluster formation. 

The authors thank Pieter van Dokkum for making the data 
available in reduced form, which has been invaluable for this 
analysis. 
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18 Weak lensing analysis ofMSl 054-03 

APPENDIX 

ERROR ANALYSIS 

The weak lensing signal is obtained by averaging the shape measurements of a number of faint galaxies. The noise in the shape 
estimate differs from object to object, and therefore it is important to weight the measurements properly. 

For bright galaxies the error on the shear measurement is dominated by the intrinsic ellipticity distribution of the galaxies. At 
fainter magnitudes, shot-noise increases the uncertainty in the measurements. As these galaxies require large seeing corrections, the 



measure ment e rror on the distortion can be quite high. This is illustrated in figure A20 



Figure A20i shows the observed scatter in the polarization as a function of apparent magnitude in the F814W band. As fainter 



objects suffer more from the circularization by the PSF one would expect the scatter to decrease towards fainter magnitud es. H owever 



the increasing contribution by shot noise (indicated by the solid line) cancels out the effect of the circularization. Figure A2C b shows 

the observed scatter in the ellipticities (shear) of the objects. The scatter at faint magnitudes increases rapidly due to the combination 

of shot noise and large PSF corrections. The solid Une indicates the expected scatter due to the combination of shot noise and PSF 

correction. 

We estimate the weighted mean of the distortion using: 

{S^) = ^^^P^, (Al) 

where the weight w„ is the inverse of the variance of the distortion: 

1 p2 

where (7q) is the scatter due to the intrinsic ellipticity of the galaxies, P^ is the pre-seeing shear polarizability and (Ae^)^^ is the 
error estimate for the polarization, which is derived below. 

This procedure, however, gives only the optimal estimate for the average ellipticity, but does not account for the fact that the 
nearby bright galaxies are lensed less efficiently than faint high redshi ft gal axies. To get the best lensing estimate it is therefore 



also necessary to weight with (3 as well. This is demonstrated in figure A21a, where we show the weight function as a function of 
apparent magnitude (solid line). The dashed line shows the weight function hen we include the strength of the lensing signal. This 
effectively decreases the weight of bright (nearby) galaxies and increases the relative weight of the faint (distant) galaxies. 



It is also interesting to see which magnitude range contributes most to the measurement of the distortion. In figure |A21 b we plot 

each bin 



the weight function for each magnitude bin times the number of galaxies in that bin. This gives the total cont ributi on of 



to the measurement of the distortion. Although the bright galaxies have the highest weight function (cf. figure A21 1), they are not 



so numerous to contribute significantly to th e lensing signal. At the faint end both incompleteness and increasing noise causes the 



contribution to decrease rapidly. Figure A21b also shows that including the strength of the lensing signal (/3) changes the profile only 



slightly. It shows that galaxies with FS14W magnitudes between 24 and 26.5 contribute most to the lensing signal. 

Estimating the measurement error on the polarization 

Here we derive an expression to estimate the error on the polarization measurement, using the observations. 
The observed quadrupole moments are given by: 

7°/^ = / flf2«,x,W(i)/obs(i), (A3) 

which are combined to form the two-component polarization. 

I\\ — 122 27i2 

ei = and ea = (A4) 

hi+hi hi+hi 

Due to noise, the true image /ti-ue(.^ is changed to the observed image fobsix): 

/obs {x) = /ti-ue (i) + ab{x) + d(x) \/ f{x) , (A5) 

where we take {b{x)^)2 = {d{xf-)'- = 1. The first term in equation A5 gives the true image, the second term is the noise in the 
background, and the thkd contribution comes from the shot noise due to the image. We therefore observe quadrupole moments given 
by: 

lf°' = Ifp + I d^xXiXjW(x) (abix) + d(x) ^/fi^) (A6) 

As {b{x)) = {d{x)) = 0, this shows that (1°^^) = /,'""'. Thus noise does not introduce a bias in the average measurement of the quadrupole 
moments. However, it does introduce a neglegible bias in the polarization. We now can estimate the variance in the observed 
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Fig. A20. — (a) Plot of the observed scatter in the polarization 
as a function of apparent magnitude in the F814W band. The 
scatter is fairly constant with apparent magnitude. The PSF 
circularizes the images of the faint galaxies, but the increas- 
ing scatter due to shot noise cancels out this effect, (b) Plot of 
the observed scatter in the ellipticities (shear) of the galaxies. 
At bright magnitudes the scatter is dominated by the intrinsic 
ellipticities of the galaxies. At fainter magnitudes the scatter 
increases rapidly due to increasing shot noise and large PSF 
corrections. The contribution of the shot noise to the scatter is 
indicated by the solid line. 

quadrupole moments using equation A6 and noting that {b{x)b(y)) ■ 
{b(x)b(y)6ix-y)), {d(x)d(y)) = {d(x)d{y)6{x-y)), and {b{x)d(y)) = 
0, where 6{x-y) is the Dirac delta function. This yields; 



(AlijAhi) = / d^xxiXjXkX,W^(x) {a^+f(x)) (A7) 
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Fig. A21. — (a) Weight function as a function of apparent 
F814W magnitude (solid line). For bright galaxies the weight 
function is dominated by the intrinsic ellipticities of the galax- 
ies. At fainter magnitudes, the weight function decreases 
rapidly due to the increasing measurement error on the distor- 
tion. Including the strength of the lensing signal (/3) yields a dif- 
ferent weight function (dashed line), which lowers the weight of 
bright (nearby) galaxies and thus increases the relative weight 
of faint (distant) galaxies, (b) Plot of the weight function times 
the number of objects in each magnitude bin as a function of ap- 
parent F814W magnitude (solid line). This plot shows which 
magnitude bin contributes most to the signa l. Al though bright 
galaxies have the highest weight (cf. figure i\21[ ), they are not 
so numerous that they contribute significantly to the lensing sig- 
nal. Including the strength of the lensing signal (fi) changes the 
profile mainly at the bright end, but the differences are small. 
The lensing signal is dominated by galaxies with F814W mag- 
nitudes between 24 and 26.5. 

We now want to estimate the errors on the measurements of 
the polarizations e\ and e2 and express them in terms of the 
errors on the quadrupole moments. This finally gives: 
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Weak lensing analysis of MS 1054-03 
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Fig. A22. — Left panel: ElUpticity introduced by constructing an interlaced images from images with imperfect offsets. The reference 
star (a Tiny Tim PSF) is centered exactly on a pixel. This image is shifted (Ajc, Ay) and an interlaced image is created, assuming 
the shift was (0.5, 0.5) pixel; Right panel: The same procedure was repeated for the case in which the reference star is placed at 
the intersection of four pixels. The sticks indicate the direction and the size of the introduced ellipticity. If the shifts are half pixel, 
the measured ellipticity is correct. The offsets of our observations for F814W (open circles) and F606W (solid triangles) are also 
indicated. For one F814W pointing the offset is such that results of the interlaced images are biased. 
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From simulations and real images we find that the errors on 
both e\ and 62 are similar. This estimate for the error on the 
polarization can be calculated from the observations and allows 
us to properly weight the signals of the galaxies using equa- 
tion A2. 

THE EFFECTS OF INTERLACING 

The images of small, faint galaxies in WFPC2 observations 
are poorly sampled. HFKS98 showed that as a result the small- 
est objects cannot be corrected reliably for the PSF effects. 

However, several techniques can be used to improve the 
sampling of WFPC2 observations by dithering the exposures 
(Fruchter & Hook 1998). When the shifts are half pixel in both 
X and y (modulo integer pixel shifts) one can construct an in- 
terlaced image (Fruchter & Hook 1998) which has a V2 better 
sampling. 

Interlaced images were constructed from the observations of 
MS 1054-03. Doing so, we assumed that the offsets were ex- 



actly half pixel. Due to pointing errors and the camera distor- 
tion the real offsets are slightly different. 

In this section we examine the effect of imperfect offsets on 
the shape measurements. To do so, we examine the change in 
the shape of a model PSF. Although interlacing is not a convo- 
lution, it affects large objects (as these are well sampled even 
before interlacing) less than small objects. By examining the 
effects on the PSF, we study a worst case scenario. The effects 
should be less for most of the galaxies used in the weak lensing 
analysis. 

We use a ten times oversampled model PSF calculated using 
Tiny Tim. This PSF is shifted and rebinned to WFPC2 sampling 
and convolved with the pixel scattering function, in order to 
realistically mimic the WFPC2 PSF. 

We start with a reference image in which the star is centered 
exactly on a pixel. This image is shifted (Ax, Ay) and an in- 
terlaced image is created assuming the shift was exactly (0.5, 
0.5) pixels. The resulting image is analysed and the measured 
ellipticities are compared to the correct value. 

The change in el liptic ity as a function of the applied offset is 
presented in figure A22 (left panel). When the reference star is 
placed at the intersectio n of f our pixels the results are different, 
as is indicated in figure A.22 (right panel). These two extreme 
cases shoul d cov er the full range of possible configurations. 

In figure \22 we also plot the observed offsets for our data. 
The open circles represent the offsets of the F814W images, 
and the solid triangles show the offsets for the F606W images. 
Only one F814W pointing has such a deviating offset that the 
shape measurements cannot be used for our weak lensing anal- 
ysis. In all other images the offsets are such that no significant 
bias is introduced by the interlacing. 
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Fig. B23. — Ellipticity introduced by imperfect offsets and WFPC2 camera distortion, (a): The introduced ellipticities for the worst 
offset in the F814W filter (Ajc, Ay) = (5.71,5.89) as a function of position. In this case it is clear that interlacing introduces a bias in 
the observed shapes, and therefore these data (the F814W image of the first pointing) are not used in the weak lensing analysis, (b): 
The introduced ellipticity for the second worst offset (5.53, 5.60). In this case the bias is very small. Due to the telescope distortion 
the effect is somewhat larger at the edges of the chips. The sticks indicate the direction and the size of the introduced ellipticity. 



Even if the offset is perfect in the centre of the image, tele- 
scope distortions will result in imperfect offsets at the edges. 
We simulated the amplitude of the effect using the coefficients 
for the telescope distortions from Holtzman et al. (1995). 

Figure |B23| shows the introduced ellipticity for our worst off- 
set (5.71, 5.89). The bias in the observed ellipticities is clear in 
this case. These data (the F814W image of the first pointing) 
were not used in the weak lensing analysis. The second worse 
offset (5.53, 5.60) is also shown in this figure. The ellipticity 
introduced in this case is small. Only at the edges of the chips 
the camera distortion introduces a small effect. As the effect is 
even smaller for galaxies, we feel certain that we can use the 
interlaced images for our weak lensing analysis. 



